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We review supersymmetric models where R-parity is broken either ex- 
plicitly or spontaneously. The simplest unified extension of the MSSM with 
explicit bilinear R-Parity violation provides a predictive scheme for neu- 
trino masses and mixings which can account for the observed atmospheric 
and solar neutrino anomalies. Despite the smallness of neutrino masses 
R-parity violation is observable at present and future high-energy collid- 
ers, providing an unambiguous cross-check of the model. This model can 
be shown to be an effective model for the, more theoretically satisfying, 
spontaneous broken theory. The main difference in this last case is the 
appearance of a massless particle, the majoron, that can modify the decay 
modes of the Higgs boson, making it decay invisibly most of the time. 

PACS numbers: 14.60.Pq, 11.30.Pb, 12.60.Jv 

1. Introduction 

Despite the tremendous effort that has led to the discovery of neutrino 
mass III, |j, y| the mechanism of neutrino mass generation will remain open 
for years to come (a detailed analysis of the three-neutrino oscillation pa- 
rameters can be found in jj|). The most popular mechanism to generate 
neutrino masses is the seesaw mechanism [!, H,ui H, HI • Although the seesaw 
fits naturally in SO(10) unification models, we currently have no clear hints 
that uniquely point towards any unification scheme. 

Therefore it may well be that neutrino masses arise from physics having 
nothing to do with unification, such as certain seesaw variants [T^j, and 
models with radiative generation jlH ll^ . Here we focus on the specific 
case of low-energy supersymmetry with violation of R-parity, as the origin 
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of neutrino mass. R-parity is defined as R p = (_1) 3B + L + 25 with S, B, L 
denoting spin, baryon and lepton numbers, respectively [Hi ]. 

In these models R-parity can be broken either explicitly or sponta- 
neously. In the first case we consider the bilinear R-parity violation model, 
the simplest effective description of R-parity violation T4]. The model 
not onl y a ccounts for the observed pattern of neutrino masses and mix- 
ing 2) III 13- 3- but also makes predictions for the decay branching ra- 



tios of the lightest super symmetric particlejl9l. 20l . 21 , z| from the current 
measurements of neutrino mixing angles 4]. In the second case R-parity 
violation takes place "a la Higgs", i.e., spontaneously, due to non-zero sneu- 
trino vacuum expectation values (vevs) [23l 24|, 2f|. In this case one of the 
neutral CP-odd scalars is identified with the majoron, J. In contrast with 
the seesaw majoron, ours is characterized by a small scale (TeV-like) and 
carries only one unit of lepton number. In previous studies |2fJ, |27l, 2g, 2{| 
it was noted that the spontaneously broken R-parity (SBRp) model leads 
to the possibility of invisibly decaying Higgs bosons, provided there is an 
SU(2) (g) U(l) singlet superfield <I> coupling to the electroweak doublet Hig- 
gses, the same that appears in the NMSSM. We have reanalyzed [3(| this issue 
taking into account the small masses indicated by current neutrino oscilla- 
tion data 0]. We have shown explicitly that the invisible Higgs boson decay 
Eq. CO, 

h-> J J (1) 

can be the most important mode of Higgs boson decay. This is remark- 
able, given the smallness of neutrino masses required to fit current neutrino 
oscillation data. 



2. Bilinear R-Parity Violation 

2.1. The Model 

Since BRpV has been discussed in the literature several times 0, h], 3, 
3~D . 32l . 33l ] we will repeat only the main features of the model here. We will 
follow the notation of jltl IriT ] . The simplest bilinear Ijt p model (we call it 
the I/l p MSSM) is characterized by the superpotential 

W = Wmssm + W 9p (2) 
In this equation WmSSM is the ordinary superpotential of the MSSM, 



W = e ab 



hgQfUjH* + h%Q\D 3 H a d + 4^RjH a d - ixH a d H b u (3) 



where i,j = 1,2,3 are generation indices, a, b = 1,2 are SU(2) indices. We 
have three additional terms that break R-parity, 

W^ p = e ab e^H b u . (4) 
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These bilinear terms, together with the corresponding terms in the soft 
supersymmetric (SUSY) breaking part of the Lagrangian, 



Ksoft = K)oft SM + e abBi€iL°;Hu 



(5) 



define the minimal model, which we will adopt throughout this paper. The 
appearance of the lepton number violating terms in Eq. @ leads, in gen- 
eral, to non-zero vacuum expectation values for the scalar neutrinos (z/j), 
called Vi in the rest of this paper, in addition to the VEVs v u and Vd of 
the MSSM Higgs fields H® and H^. Together with the bilinear parameters 
£j, the Vi induce mixing between various particles which in the MSSM are 
distinguished (only) by lepton number (or R-parity). Mixing between the 
neutrinos and the neutralinos of the MSSM generates a non-zero mass for one 
specific linear superposition of the three neutrino flavor states of the model 
at tree-level while 1-loop corrections provide mass for the remaining two 
neutrino states 0, 0] . 

2.2. Tree Level Neutral Fermion Mass Matrix 

In the basis V> 0T = (— iA', — Hj, v e , v T ) the neutral fermions 
mass terms in the Lagrangian are given by 



C m = -^(ij ) T M N yj + h.c. 



(6) 



where the neutralino/neutrino mass matrix is 

M N = 



M x o m T 




in 



(7) 



with 



Mi 




y 

L 29 >v u 



M 2 



-29 v d 2~gVd 



-\g'vd 
\gvd 
o 



\g'v u 
-\gv u 



m 



Ol 
«3 



(8) 



where a% = (— ^g'vi, ^gVi, 0, ej). This neutralino/neutrino mass matrix is 
diagonalized by 

M*M N N~ l 



diag(m x o , m x o , m x o , m x o , m ui , m V2 , m Vs ) 



(9) 
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2.2.1. Approximate Diagonalization at Tree Level 

If the Ijlp parameters are small (we will show below that this is indeed 
the case), then we can block-diagonalize Mn approximately to the form 
diag(m cff ,.M x o) 



m cS = -m ■ M x lm T 



with 



Miff 2 + M 2 g 
4 det(.MV 



/2 



( kl 



AeA,, 
A 2 



\ A e A T A^A T 



A e A r 
A M A r 
A 2 



Aj = (IVi + V d €i . 

The matrices N and V u diagonalize M x o and m e s 

N*M x oN< = diag(m x p) ; Vj m e f[V v = diag(0, 0, m v 

where 



Tr(m eS ) 



Miff 2 + M 2 ff /2 ^ |2 



(10) 



(11) 



(12) 



(13) 



4 det(M x o) 

So we get at tree level only one massive neutrino, the two other eigen- 
states remaining massless. The tree level value will give the atmospheric 
mass scale, while the other states will get mass at one loop level. Therefore 
the BRpV model produces a hierarchical mass spectrum. We will show below 
how the one loop masses are generated. 

2.3. One-Loop Neutrino Masses and Mixings 
2.3.1. Definition 

The self-energy for the neutralino/neutrino is [lii ]. 



{j> [Pl4 + PRsg] - [p L uf- + p R n%] } (14) 




Then the pole mass is, 



M pole = M DR (/i/?)+AM) 



1.7 



(15) 



with 



AM,', 



\ ( n S ( m l ) + n S (m ? ) ) - \ (jn x o Sg (m • ) + m x o Eg (rrijfj 



A=0 



(16) 
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where 

E^ = |(E i + E il ) ; n v = i(n L + n fl ) (17) 

an the parameter that appears in dimensional reduction is, 

2 

A = — - - 7i; + ln47r . (18) 

Explicit expressions can be found in [lsl ]. 
2.3.2. One—loop results for the masses 




Fig. 1. a) One- loop neutrino masses, b) Atmospheric angle as a function of 
lAJ/v'Af+A?. 

The BRpV model produces a hierarchical mass spectrum for almost all 
choices of parameters. The largest mass can be estimated by the tree level 
value using Eq. (|13|). Correct Am 2 TM can be easily obtained by an appro- 
priate choice of |A|. The mass scale for the solar neutrinos is generated 
at l-loop level and therefore depends in a complicated way in the model 
parameters. We will see below how to get an approximate formula for the 
solar mass valid for most cases of interest. Here we just present in Fig. ^ 
a), for illustration purposes, the plot of the three eigenstates as a func- 
tion of the parameter e 2 /|A|, for particular values of the SUSY parameters, 
itiq = \x = 500 GeV, tan/3 = 5, B = —A = tuq. For the R-parity parameters 
we took |A| = 0.16 GeV, 10 * A e = A M = A T and €\ = e 2 = e 3 . 

2.3.3. The mixings 

Now we turn to the discussion of the mixing angles. As can be seen from 
Fig. |I]a), if e 2 /|A| <C 100, then the l-loop corrections are not larger than 



G 
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the tree level results. In this case the flavor composition of the 3rd mass 
eigenstate is approximately given by 

U a3 « A Q /|A| (19) 

As the atmospheric and reactor neutrino data tell us that Un — > v T oscilla- 
tions are preferred over — ► v e , we conclude that 

A e < A M ~ A r (20) 

are required for BRpV to fit the data. We cannot set all the A, equal, because 
in this case U e3 would be too large contradicting the CH00Z, as shown in 
Fig. da). 
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Fig. 2. a) Atmospheric angle as a function of |A M |/ ^/A^ + A£. b) U^ 3 as a function 



Of lAel/y/AJ+A?. 



We have shown that even a very small deviation from universality (less 
then 1%) of the soft parameters at the GUT scale allows for 



e e , A e 
A 



(21) 



A' 



Then we can have at the same time small U^ 3 determined by A e /A^ as in 
Fig. [2]a) and large tan 2 (# SOL ) determined by e e /e M as in Fig. [21 b). After 
the Kamland and SM0 salt results, this is the only scenario consistent with 
the data. 



2-4- Approximate formulas for the solar mass and mixing 

In all the previous analysis we used a numerical program to evaluate 
the one- loop masses and mixings. It is however desirable to have analytical 



GustavoFest-Corfu2005 printed on February 7, 2008 



7 



approximate results that can give us quickly the most important contribu- 
tions. We have identified that these are the bottom-sbottom loops and the 
charged scalar-charged fermion loops. Then, by expanding in powers of the 
small R-Parity breaking parameters 6j, we get approximate formulas for the 
solar mass scale and mixing angle as we will explain below. 



2.4.1. Bottom-sbottom loops 

As an example, we give the result from the bottom-sbottom loop [l5j |. 



AM 



j—^2s- b c- b h b AB 



+ a 3 6(ei<5j 3 + €j5 i3 ) |A| 



2 . a L aB.\ x 



(22) 



where are the in the basis where the tree level neutrino mass matrix is 
diagonal, 



v mf e 



.i ■ 



(23) 



the a,i are functions of the SUSY parameters, and 



AB^ = B (0,m 2 b ,ml) ~ fi o(°- m '' m l) ■ 



(24) 



where the Bq are the Passarino-Veltman loop functions |2J]. The different 
contributions can be understood as coming from different types of insertions 
as shown in Fig. El In this figure open circles correspond to small R-parity vi- 
olating projections, full circles to R-parity conserving projections, and open 
circles with a cross inside to mass insertions which flip chirality. With this 
understanding one can make a one to one correspondence between Eq. (|22j) 
and Fig. |21 




"J W,B i , , B „ v * 
■i y y 




-o- 



a s |A|* j3 +&ej b a s \A\S is + bii (a 2 , ai)|A|5 j3 b ai|A|5 i3 
Fig. 3. Bottom-sbottom loop and different types of insertions. 
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2.4.2. Simplified approximation formulas 

In the basis where the tree-level neutrino mass matrix is diagonal the 
mass matrix at one-loop level can be written as 

/ 



V^ T m u V^ 



where 



co 



cieiei ciei£2 

Ci? 2 ?l Ci? 2 ?2 
y Cie 3 ?i Ci? 3 ?2 



M l9 2 + M 2 g u 
4det(A4 x o) 



cieie 3 



\ 



c |A| 2 



+ • 



+ cie 3 e 3 J 



(25) 



Cl = T^ s ' m{26 ~ b)mbABb ° 2bl 7? 



(26) 



The dots in Eq. (|25|) correspond to the terms that are not proportional to 
the ii x ij structure, as can be seen from Eq. (|22|) . Assuming that the 
bottom-sbottom loop dominates, the e« x ij structure is dominant, and the 
matrix can be diagonalized approximately under the condition 



Q)|A| 2 



Then we get 



and 



16vr 2 



sin(2fc)m 6 AB l 







Il + ii 

M 2 



2 

tan SOL = 

4 



(27) 



(28) 



(29) 



The results for the masses are presented in Fig. |1J On the left panel we 
have the data set where the neutralino is the LSP while on the right panel 
the LSP is the stau. We see that the agreement is fairly good, particularly 
in the region allowed by the present data. 



2.5. Probing Neutrino Mixing via SUSY Decays 

After having shown, in the previous sections, that the BRpV model pro- 
duces an hierarchical mass spectrum for the neutrinos that can accommo- 
date the present data for neutrino masses and mixings, we now turn to ac- 
celerator physics and will show how the neutrino properties can be probed 
by looking at the decays of super symmetric particles. 
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Fig. 4. Comparison between the simplified analytical formulas and numerical results 
for the solar mass. In the left panel the LSP=x° while in the right panel the LSP=f . 
The red (dark) band corresponds to the latest neutrino data 



If R-parity is broken the LSP will decay. If the LSP decays then cosmo- 
logical and astro-physical constraints on its nature no longer apply. Thus, 
within R-parity violating SUSY, a priori any superparticle could be the LSP. 
In the constrained version of the MSSM (mSUGRA boundary conditions) one 
finds only two candidates for the LSP, namely the lightest neutralino or one 
of the right sleptons, in particular the right scalar tau. Therefore we will 
consider below these two cases. However, if we depart from the mSUGRA 
scenario, it has been shown recently p3] that the LSP can be of any other 
type, like a squark, gluino, chargino or even a scalar neutrino. 



2.5.1. Probing Neutrino Mixing via Neutralino Decays 

If the neutralino is the LSP and R-parity is broken, it will be unstable 
and it will decay. It was shown 1 in Ref . [2(| , that the neutralino decays well 
inside the detectors and that the visible decay channels are quite large. We 
have shown in the previous sections that the ratios |Aj/Aj| and j / | were 
very important in the choice of solutions for the neutrino mixing angles. 
What is exciting now, is that these ratios can be measured in accelerator 
experiments. In the left panel of Fig. \5\ we show the ratio of branching 
ratios for semileptonic neutralino decays into muons and taus: BR(x 
Hq'q)/BR(x — > T~q'q) as function of tan 2 # ATM . We can see that there is a 
strong correlation. The spread in this figure can in fact be explained by 
the fact that we do not know the SUSY parameters and are scanning over 



The relation of the neutrino parameters to the decays of the neutralino has also been 
considered in Ref. [35||. 
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the allowed parameter space. This is illustrated in the right panel where we 
considered that SUSY was already discovered with the following values for 
the parameters, 

M 2 = 120 GeV, n = m = 500 GeV, tan (3 = 5, A = -500 GeV. (30) 




left for random SUSY values and on the right for the SUSY point of Eq. H30|) . 



2.5.2. Probing Neutrino Mixing via Charged Lepton Decays 

After considering the case of the LSP being the neutralino, for complete- 
ness, we have also studied jis| the case where a charged scalar lepton, most 
probably the scalar tau, is the LSP. We have considered the production and 
decays of f, e and /I, and have shown that also for the case of charged 
sleptons as LSPs they will decay well inside the detector. We can correlate 
branching ratios with neutrino properties. This is shown in Fig. where the 
branching ratios of the scalar tau show a strong correlation with the solar 
angle. 

2.5.3. Other LSP decays 

As we discussed in the introduction to this section, if we depart from 
the mSUGRA scenario, then the LSP can be of other type, like a squark, 
gluino, chargino or even a scalar neutrino, as it has been shown recently in 
Ref. [22|. As the decays of these LSP's will always depend on the parameters 
that violate R-parity and induce neutrino masses and mixings, it is possible 
to correlate branching ratios to the neutrino properties. This is shown in 
Fig. Qtaken from Ref 22] . On the left panel we consider the case of chargino 
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tan 2 Q tan 2 

Fig. 6. Ratios of branching ratios for scalar tau decays versus tan 2 # Q . The left 
panel shown all data points, the right one refers only to data points with €2/63 
restricted to the range [0.9,1.1]. 




(A2/A3) 2 tan 2 Sol 
Fig. 7. Left panel: ratios of branching ratios for chargino decays versus (A2/A3) 2 ; 
Right panel: ratios of branching ratios for squark decays versus tan 2 6*SOL. These 
plots were taken from Ref. [22| . 

decays and show BR(x + — > /j,cc)/BR(x + — ► tcc) as a function of the ratio 
(A2/A3) 2 . As we saw, this last ratio is correlated with the atmospheric 
angle. So, by looking at the chargino decays one can test the atmospheric 
mixing angle. On the right panel we show a very strong correlation for the 
solar mixing angle obtained [22| with the decays of squarks. Many other 
correlations can be obtained making the model over constrained. So, in 
summary, no matter what supersymmetric particle is the LSP, measurements 
of branching ratios at future accelerators will provide a definite test of the 
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BRpV model as a viable model for explaining the neutrino properties. 

3. Spontaneous Broken R-Parity Model 

3.1. The Model 

The most general superpotential for the SBRp model is [2^l . 24l . 2fil . H3| 

W = e ab (4QtUjH b n + h%Q\D 3 H a d + h%L%R a d 
+h^LfV]H b u - mH h u - (h H2H b u + 5 2 )$ 



+ tf'*^ + M^S, + -M $ $ 2 + ^ (31) 



The first three terms together with the fi term define the R-parity conserving 
MSSM, the terms in the last row only involve the SU(2) <g) U(l) singlet super- 
fields (v c i, Si,<&) carrying a conserved lepton number assigned as (—1, 1,0), 
while the remaining terms couple the singlets to the MSSM fields. To this we 
will have to add the soft SUSY breaking terms in the usual way [30j. 

3. 2. Pattern of spontaneous symmetry breaking 

The spontaneous breaking of Rp is driven by nonzero vevs for the scalar 
neutrinos. The scale characterizing Rp breaking is set by the isosinglet 
vevs, 

(32) 

We also have very small left-handed sneutrino vacuum expectation values 

= ^ (33) 
The electroweak breaking is driven by 

{Hu) = W {Hd) = T2 (34) 

with v 2 = v 2 + v\ + Yli v Li an d m w = ^T~- The spontaneous breaking 
of R-parity also entails the spontaneous violation of total lepton number. 
This implies that the majoron 



J = Im 



(v u H u - v d H d ) + 2^ -TfVi + ~T7 ~ ~T7 V 



(35) 



remains massless, as it is the Nambu-Goldstone boson associated to the 
breaking of lepton number. 
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3.3. Neutrino- Neutralino-Singlino mass matrix 

For simplicity we choose only one generation of singlet superfields. Then 
in the basis (— iX', —i\ 3 , Hd, H u , u e , u^, u T , v c , S, $) the neutrino/neutralino 
mass matrix is 10 x 10, 



M N 



( M x° 


m x o u 


m x o 


T 

m n 

X°v 





m D 


T 


T 
m D 























M„c S 








where M Y o is the usual MSSM neutralino mass matrix, and 



T 



( 

V 



~\g'vLe 


\gVLe 
















-\9'vLt 


\gvLr 







m D )i = 


T2 Kvu 







with an effective "bilinear R-parity" parameter 

1 



V2 



K V R 



(36) 



(37) 



(38) 



(39) 



The other block matrices in Eq. ()36|) can found explicitly in Ref.[3 

3.4- The effective neutrino mass matrix 

The exact diagonalization of Eq. ()36|) can only be done numerically. 
However due to the smallness of the neutrino masses, an effective neutrino 
mass matrix can be cast into a very simple form 



{mf u )ij = oAjAj + b(eiAj + e^Aj) + cejCj 



(40) 



where the model dependence in hidden in the parameters a, b, c. The effec- 
tive bilinear R-parity violating parameters ej and Aj, are given in Eq. (|39|) 
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and Eq. respectively. This equation resembles very closely the result 
for the BRpV model once the dominant 1-loop corrections are taken into ac- 
count. The tree-level result of the explicit bilinear model can be recovered 
in the limit Mr, M<j> — > oo. 

g 2 M 1+ g' 2 M 2 

a = — , 6 = 0, c = (41) 

4DetM x o ' K ' 

In this limit only one non-zero neutrino mass remains. The relative size of 
the coefficient c compared to the corresponding 1-loop coefficient dictates 
if the 1-loop corrections or the contribution from the singlet fields are more 
important. Both extremes can be realized in our model. Large branching 
ratios of the Higgs into invisible final states require sizable values of h and 
/iq. Then the "singlino" contribution dominates. 



3.5. The Neutral Scalar Bosons Mass Matrix 

The 8x8 scalar mass matrix is a symmetric matrix that in the basis 
(H%, H®, Di, <3?, S, v c ) takes the form, 







M s2 ~ 

HL 






M s2 = 








(42) 




_ </ 


M l7 


Mf s _ 





written in terms of block matrices that can be found in Ref.|30|. The 
pseudo-scalar mass matrix has a similar form. In practice they are quite 
complicated and have to be diagonalized numerically. The most important 
feature of the pseudo-scalar mass matrix is the existence of two massless 
Goldstone bosons, corresponding to the spontaneous breaking of gauge sym- 
metry (the usual Goldstone Go "eaten" by the Z°) and to the spontaneous 
breaking of R-parity, the majoron J. In the basis 

A^ = (Hf,H° u I y I ,i> 2I ,u SI ,^,S 1 ,^) (43) 

we have 

G = (N v d ,-N v u ,N v L1 ,Nov L2 ,N v L3 ,0,0,0) 

J = (-Niv d ,N 1 v u ,N 2 v L1 ,N 2 v L 2,N 2 v L3 ,0,N 3 vs,-N 3 v R ) (44) 

where the iVj are normalization constants. It can easily be checked that 
they are orthogonal Gq ■ J = 0. 
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3.6. Higgs Boson Production 

Supersymmetric Higgs bosons can be produced at an e + e _ collider via 
the so-called Bjorken process, 



Chzz = ^(v^Gf) 1/2 M| Z^^tf; (45) 

i=l 

where the r/j are combinations of the doublet scalars vevs and rotation ma- 
trices: 

^v4 + 74+E% (^) 

In comparison with the SM the coupling of the lightest CP-even Higgs boson 
to the Z is reduced by a factor 

r) = rji < 1. (47) 
In the MSSM we have r] = sin(/3 — a). 



3. 7. Higgs Boson Decay 
We are interested here in the ratio 

_ Tjh -> J J) 
Rjb ~ T(h -+ 66) (48) 

of the invisible decay to the SM decay into b-jets. These decay widths are 
easily obtained [3oT |. 



and 

r / \ 2i 3 / 2 



T(h -» 66) = g y (-Rii) m ft m£ 



I — ■ 1 ' ' 2 



(50) 



We want to look at situations where Rjb > 1, but that at the same time the 
Higgs boson is produced at a reasonable rate, which means that r\ parameter 
should not be too small. 
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3.8. Scanning Strategy 

As the model has a large number of parameters we have to follow a 
scanning strategy that can put in evidence the possibility of the Higgs boson 
decaying invisibly. This is achieved by fixing as many parameters as we can. 
Without loss of generality we then fix the MSSM parameters to the SPSla 
benchmark point: 

m = lOOGeV, m 1/2 = 250GeV, tan = 10 

Aq = -lOOGeV, \i < (51) 

For the singlet parameters we choose to start with 

V R = v s = v $ = — 150GeV 

M R = -M* = 5 = 10 3 GeV 

h = 0.8 , ho = -0.15 , A = 0.1 (52) 

and then perform variations around these values in a controlled way. Finally 
the explicit bilinear parameters, and Aj, are fixed approximately such that 
neutrino masses and mixing angles are in agreement with experimental data. 



3.9. Numerical results: General Superpotential 

In Fig. |H] we show the plot of Rjf, as a function of if' . On the left panel 
we kept h constant and varied continuously the other singlet parameters, 
while on the right panel the vevs were kept. We see that it is indeed possible 
to achieve, at the same time, both large branching ratio into the invisible 
mode (h — ► J J), and large production cross sections. In Ref. |3fl] a complete 
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Fig. 8. Rjb as a function of n 2 : a) For constant values of h. b) For constant values 
of v R = v s - 



analysis of the parameter space where this happens is presented. 
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4. Conclusions 



The BRpV model is a simple extension of the MSSM that leads to a very 
rich phenomenology. We have calculated the one-loop corrected masses and 
mixings for the neutrinos in a completely consistent way, including the RG 
equations and correctly minimizing the potential. We have shown that it 
is possible to get easily maximal mixing for the atmospheric neutrinos and 
large angle MSW, as it is preferred by the present neutrino data. We have 
also obtained approximate formulas for the solar mass and solar mixing 
angle, that we found to be very good, precisely in the region of parameters 
favored by this data. 

We emphasize that the LSP decays inside the detectors, thus leading 
to a very different phenomenology than the MSSM. No matter what super- 
symmetric particle is the LSP measurements of branching ratios at future 
accelerators will provide a definite test of the BRpV model as a viable model 
for explaining the neutrino properties. 

We have also discussed the possibility of invisibly decaying Higgs boson 
in the context of the SBRp model. One of the neutral CP-odd scalars in 
this model, is the Nambu-Goldstone boson associated to the breaking of 
lepton number, the majoron J. Neutrino masses and mixings can easily 
be explained. In contrast to the MSSM, where the Higgs boson can decay 
invisibly only to supersymmetric states, in our case the Higgs boson can 
decay mainly invisibly through the decay h — ► J, J. As a result, our analysis 
indicates that invisibly decaying Higgs bosons should be an important topic 
in the agenda of future accelerators, such as the Large Hadron Collider and 
the Next Linear Collider. 
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